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Abstract 

A new narrow 3-atom loss resonance associated with an Efimov trimer crossing the 3-atom 
threshold has recently been discovered in a many-body system of ultracold 6 Li atoms in the three 
lowest hyperfine spin states at a magnetic field near 895 G. O'Hara and coworkers have used 
measurements of the 3-body recombination rate in this region to determine the complex 3-body 
parameter associated with Efimov physics. Using this parameter as the input, we calculate the 
universal predictions for the spectrum of Efimov states and for the 3-body recombination rate 
in the universal region above 600 G where all three scattering lengths are large. We predict an 
atom-dimer loss resonance at 672 ±2 G associated with an Efimov trimer disappearing through 
an atom-dimer threshold. We also predict an interference minimum in the 3-body recombination 
rate at 759 ±1 G where the 3-spin mixture may be sufficiently stable to allow experimental study 
of the many-body system. 

PACS numbers: 31.15.-p, 34.50.-s, 67.85.Lm, 03.75.Ss 
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I. INTRODUCTION 



Ultracold atomic gases allow the experimental study of many-body physics in systems 
in which the fundamental interactions are theoretically understood and experimentally con- 
trollable. There are many possibilities for the atomic constituents of the many-body system: 
they can be identical atoms that are all in the same spin state or there can be multiple species 
of constituents that are either different spin states of the same atom, different isotopes of an 
atom, or atoms of different elements. If the atoms are identical fermions in the same spin 
state, the many-body system is essentially an ideal gas. If the atoms are fermions with two 
spin states, there is a well-behaved zero-range limit in which the interactions are completely 
determined by the scattering length a. The ground state of the many-body system with 
equal populations of the two spin states is a superfluid. As a is changed from a negative 
value to a positive value through ±00, the mechanism for superfluidity changes smoothly 
from the BCS mechanism (Cooper pairing of atoms) to the BEC mechanism (Bose-Einstein 
condensation of diatomic molecules). In the unitary limit in which the pair scattering length 
is infinitely large, the behavior of the system is constrained by scale invariance: invariance 
under rescaling space and time by arbitrary positive numbers A and A 2 , respectively. There 
have been extensive investigations, both experimental and theoretical, of the universal be- 
havior of the system in the crossover region [1]. 

Fermionic atoms with three spin states open up new possibilities for universal behavior. 
In the many-body system, there can be new superfluid phases and new mechanisms for 
superfluidity [2-6]. The universal behavior can be more complex, because there are three 
independent scattering lengths. However it is also qualitatively different because of the 
Efimov effect. This remarkable phenomenon was discovered by Vitaly Efimov in the 3-body 
problem for identical bosons [7]. In the unitary limit in which the three pair scattering 
lengths are all infinitely large, there is an infinite sequence of 3-body bound states called 
Efimov trimers with an accumulation point at the 3-atom threshold. The ratio of the binding 
energies of successive Efimov trimers is approximately 1/515. This geometric spectrum 
reflects discrete scale invariance: invariance under rescaling space and time by the same 
integer power of A ~ 22.7 and Aq ~ 515, respectively. Low-energy phenomena associated 
with discrete scale invariance are generally referred to as Efimov physics [8-10]. (See Ref. [10] 
for a summary of recent developments.) Efimov physics can also occur in other 3-body 
systems if at least two of the three pair scattering lengths are large. The discrete scaling 
factor depends on the mass ratios of the particles and their symmetries. In the case of 
fermions with 3 spin states for which the 3 pair scattering lengths are all large, the Efimov 
effect occurs with the same discrete scaling factor Ao as for identical bosons [11]. 

The first experimental evidence for Efimov physics in ultracold atoms was presented by 
Kraemer et al. [12]. Their atoms were 133 Cs atoms in the lowest hyperfine spin state, which 
are identical bosons. They observed resonant enhancement of the loss of atoms from 3-body 
recombination that can be attributed to an Efimov trimer crossing the 3-atom threshold. The 
dependence of the three-body recombination rate on the scattering length was first calculated 
in Refs. [13-15]. In particular, Esry, Greene and Burke predicted a resonant enhancement 
of the three-body recombination rate at negative values of the scattering length [14]. The 
universal line shape for the resonance as a function of the scattering length was first derived 
by Braaten and Hammer [16]. Kraemer et al. also observed a minimum in the 3-body 
recombination rate that can be interpreted as an interference effect associated with Efimov 
physics. In a subsequent experiment with a mixture of 133 Cs atoms and dimers, Knoop et al. 
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observed a resonant enhancement in the loss of atoms and dimers [17]. This loss feature can 
be explained by an Efimov trimer crossing the atom-dimer threshold [18]. The most exciting 
recent developments in the Efimov physics of 133 Cs atoms involve universal tetramer states. 
Platter and Hammer discovered that there is a pair of universal tetramer states associated 
with every Efimov trimer and they calculated their binding energies for limited ranges of 
the scattering length [19]. Von Stecher, D'Incao and Greene mapped out their spectrum for 
all scattering lengths and pointed out that resonant enhancement of 4-body recombination 
would provide a signature for these tetramers [20]. Ferlaino et al. recently observed the loss 
resonances from both tetramers in an ultracold gas of 133 Cs atoms [21]. 

Recent experiments with other bosonic atoms have also provided evidence of Efimov 
physics. Zaccanti et al. measured the 3-body recombination rate and the atom-dimer loss 
rate in a ultracold gas of 39 K atoms [22]. They observed two atom-dimer loss resonances 
and two minima in the 3-body recombination rate at large positive values of the scattering 
length. The positions of the loss features are consistent with the universal predictions with 
discrete scaling factor 22.7. They also observed loss features at large negative scattering 
lengths. Barontini et al. obtained the first evidence of the Efimov effect in a heteronuclear 
mixture of 41 K and 87 Rb atoms [23]. They observed 3-atom loss resonances at large negative 
scattering lengths in both the K-Rb-Rb and K-K-Rb channels, for which the discrete scaling 
factors are 131 and 3.51 x 10 5 , respectively. Gross et al. measured the 3-body recombination 
rate in an ultracold system of 7 Li atoms [24]. They observed a 3-atom loss resonance at a 
large negative scattering length and a 3-body recombination minimum at a large positive 
scattering length. The positions of the loss features, which are in the same universal region 
on different sides of a Feshbach resonance, are consistent with the universal predictions with 
discrete scaling factor 22.7. 

A promising fermionic atom in which to observe Efimov physics is 6 Li. For the three lowest 
hyperfine states of 6 Li atoms, the three pair scattering lengths approach a common large 
negative value at large magnetic fields and all three have nearby Feshbach resonances at lower 
fields that can be used to vary the scattering lengths [25]. The first experimental studies 
of many-body systems of 6 Li atoms in the three lowest hyperfine states have recently been 
carried out by Ottenstein et al. [26] and by Huckans et al. [27]. Their measurements of the 
3-body recombination rate revealed a narrow loss feature and a broad loss feature in a region 
of low magnetic field. Theoretical calculations of the 3-body recombination rate supported 
the interpretation that the narrow loss feature arises from an Efimov trimer crossing the 3- 
atom threshold [28-30] . Very recently, another narrow loss feature was discovered in a much 
higher region of the magnetic field by Williams et al. [31] and by Jochim and coworkers 
[32]. Williams et al. used measurements of the 3-body recombination rate in this region to 
determine the complex 3-body parameter that governs Efimov physics in this system. This 
parameter, together with the three scattering lengths as functions of the magnetic field, 
determine the universal predictions for 6 Li atoms in this region of the magnetic field. 

In this paper, we calculate universal predictions for various aspects of Efimov physics for 
the three lowest hyperfine spin states of 6 Li atoms. In Section II, we explain how universal 
predictions for 3-body observables can be calculated efficiently by solving coupled sets of in- 
tegral equations. We apply these methods specifically to the 3-body recombination rate and 
to the binding energies and widths of Efimov trimers. In Section III, we summarize previous 
experimental and theoretical work on 6 Li atoms in the universal region at low magnetic fields 
where all 3 scattering lengths are negative and relatively large. In Section IV, we use the 
complex 3-body parameter determined by Williams et al. to calculate universal predictions 
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for the binding energies and widths of Efimov trimers and for the three-body recombination 
rate in regions where one or more of the scattering lengths are large and positive. We pre- 
dict an atom-dimer resonance at 672 ± 2 G where an Efimov trimer disappears through an 
atom-dimer threshold. We predict an interference minimum in the 3-body recombination 
rate at 759 ± 1 G where the 3-spin mixture may be sufficiently stable to allow experimental 
study of the many-body system. We also discuss the implications of our predictions for the 
many-body physics of 6 Li atoms. We summarize our results in Section V. 

II. THEORETICAL FORMALISM 
A. Fermions with Three Spin States 

We consider a fermionic atom of mass m with three spin states that we label as types 
1, 2, and 3. We denote the scattering length of the pair ij by either = dji or a^, where 
(ijk) is a permutation of (123). We assume that the three pair scattering lengths ayi = 0,3, 
^23 = 0,1, and ai3 = a 2 are all much larger than the range of interactions. The two-body 
physics for fermions with two distinct spin states that have a large pair scattering length 
Ojj is very simple in the zero-range limit. The scattering amplitude for the pair ij with 
relative wavenumber q is (— l/a^ — iq)^ 1 . If is positive, there is a weakly-bound diatomic 
molecule with constituents i and j and with binding energy h 2 /(ma 2 j). We will refer to it as 
the (ij) dimer. We refer to the (12), (23), and (13) dimers collectively as shallow dimers. For 
realistic interactions with finite range R, there can also be deeply-bound diatomic molecules 
whose binding energies are comparable to or larger than the energy scale h 2 /(mR 2 ) set by 
the range. We refer to them as deep dimers. Their binding energies are insensitive to changes 
in the large scattering lengths. 

For the 3-atom problem, we take the zero of energy to be the scattering threshold for 
the three atoms. If is positive, the scattering threshold for an atom of type i and a (jk) 
dimer is —h 2 /(ma 2 k ). We will refer to this scattering threshold as the i + (jk) atom-dimer 
threshold. We will refer to triatomic molecules as trimers. A trimer whose constituents are 
atoms of types 1, 2, and 3 must have energy below the 3-atom threshold and below the 
% + (jk) atom-dimer threshold if ajk > 0. 

The three-body physics for fermions with three distinct spin states is more complicated 
than for two distinct spin states, because now there are three scattering lengths a\2, 023, and 
013. Depending on the signs of these scattering lengths, there can be 0, 1, 2, or 3 types of 
shallow dimers. The three-body physics is also much more intricate due to the Efimov effect. 
The unitary limit in which all three scattering lengths are infinitely large is characterized by 
discrete scale invariance with discrete scaling factor A = e 7 ^ ps 22.7, where s ~ 1-00624 
is a transcendental number. This requires a 3-body parameter that provides another length 
scale in addition to the three scattering lengths. If there were no such parameter, the system 
would have continuous scale invariance in the unitary limit. 

One simple choice for the 3-body parameter is a wavenumber /t* defined by the spectrum 
of 3-body bound states (Efimov trimers) in the unitary limit in which all three scattering 
lengths are infinite [8] . Their binding energies relative to the 3-atom threshold are 

*2 2 

Erp = A V (a 12 = a 23 = a 31 = ±00) , (1) 

m 

where is the binding wavenumber of the trimer labeled n* and n is an integer. There 
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are infinitely many of these Efimov trimers with an accumulation point at the threshold 
for three atoms of types 1, 2, and 3. The advantage of using the Efimov parameter as 
the 3-body parameter is that physical observables must be log-periodic functions of «*. In 
the zero-range limit, the Efimov spectrum also extends infinitely deep to — oo. For finite 
range R, the deepest Efimov trimer may have a binding energy E? as large as h 2 /{mR 2 ). 
However, universal predictions will be accurate only if its binding energy is significantly 
smaller, because range corrections are expected to be of order R{mE T /h 2 ) 1 / 2 . 

In the unitary limit, the Efimov states are sharp states with the spectrum in Eq. (1) only if 
there are no deep dimers in any of the three 2-body channels. If there are deep dimers, their 
inclusive effects can be taken into account by analytically continuing the Efimov parameter 
/t* to a complex value that is conveniently expressed in the form exp(ii]*/so), where 
and 77* are positive real parameters [16]. Making the substitution k* — > K*exp(ir)*/s ) on 
the right side of Eq. (1), we find that the Efimov trimers acquire nonzero decay widths \ 
The binding energies and widths are 

Erp = A V (ai2 = a 23 = a 31 = ±00) , (2) 

m 

r (n) ,2(n»-r i )2^ 2 ^sin(2^/s ) 

r T - a . (3) 

B. Three-body Recombination 

Three-body recombination is a three-atom collision in which two of the atoms form a 
dimer. In the case of three fermions in the same spin state, 3-body recombination is strongly 
suppressed at low temperature, because each pair of atoms has only P-wave interactions. In 
the case of two fermions of type % and a third atom of a different type j, two of the pairs 
have S-wave interactions with scattering length a^. The rate for 3-body recombination in the 
zero-range limit still decreases to as the energy E of the atoms approaches the threshold, 
decreasing like E if > [33] and like E 3 if < [34]. In the case of three distinct 
spin states, all three pairs of atoms can have S-wave interactions. There is no threshold 
suppression of three-body recombination if at least two of the three scattering lengths are 
large. If is large and positive, one of the recombination channels is into the (ij) dimer and 
a recoiling atom with complimentary spin k. If there are deep dimers in any of the three 
2-body channels, they provide additional recombination channels. If all three scattering 
lengths are negative, the only recombination channels are into deep dimers. 

The rate equations for the number densities rii of atoms in the three spin states are 

^Tli = -#3711712713. (4) 

The event rate constant K 3 can be separated into the inclusive rate constant K^ eop for 
recombination into deep dimers and the exclusive rate constants for recombination 

into each of the three possible shallow dimers: 

K 3 = X 3 deep + Kg" ™, (5) 

^shallow = R (12) + R (23) + ^13) (g) 

The term is nonzero only if > 0. 
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In the low-temperature limit, the rate constant K 3 and the exclusive rate constants 
can be expressed in terms of T-matrix elements for processes in which the initial state 
consists of three atoms in the spin states 1, 2, and 3 with momentum 0. By the optical 
theorem, K 3 is twice the imaginary part of the forward T-matrix element for 3-atom elastic 
scattering in the limit where the momenta of the atoms all go to 0: 



K 3 = 2ImT(0, 0, 0; 0, 0, 0). 



(7) 



The T-matrix element is singular as all the momenta go to zero, but its imaginary part is 
not. If ciij > 0, the rate constant for recombination into the (ij) dimer is the square 

of the T-matrix element for three atoms with momentum to scatter into the dimer and a 
recoiling atom multiplied by the atom-dimer phase space: 



4m 



3V3irhak 



|T fc (0,0,0;p, -p) 



\p\=2h/(V3a k ) 



(8) 



The dimer and the recoiling atom with complementary spin k both have momentum 
2h/(V3a k ). For convenience, we will switch to wavenumber variables in the remainder 
of the paper. 

The T-matrix elements in Eqs. (7) and (8) can be expressed in terms of amplitudes 
Aij(p, q; E) for the transition from an atom of type % and a complimentary diatom pair into 
an atom of type j and a complimentary diatom pair, with the two diatom pairs being the 
first to interact and the last to interact, respectively. The projection onto S-waves reduces 
the amplitude to a function of three variables: the relative wavenumber p of the incoming 
atom and diatom, the relative wavenumber q of the outgoing atom and diatom, and the 
total energy E of either the incoming atom and diatom or the outgoing atom and diatom. 
The rate constant K 3 in Eq. (7) can be expressed as 



K z = ^ aidjlmAijiO, 0; 0), 



m 



(9) 



where the sums are over i,j = 1,2,3. The exclusive rate constant in Eq. (8) for 3-body 
recombination into the (ij) dimer can be expressed as 



K. 



(ij) 



512n 2 h 



^a^ fe (0,2/(V3a fe );0) 



(10) 



where k is the complimentary spin to ij and the sum is over I — 1, 2, 3. 



C. STM Equations 

The 9 amplitudes Aij(p,q; E) satisfy coupled integral equations in the variable q that 
are generalizations of the Skorniakov-Ter-Martirosian (STM) equation [35]. To determine 
the rate constants for 3-body recombination in Eqs. (9) and (10), it is sufficient to set the 
relative wavenumber in the initial state to and the total energy to 0. The 9 coupled STM 
equations for Aij(0,p; 0) are [28] 

Ai(0,p;0) = ^ + -Y7l-<^) / dqQ(p,q;0)D k (3q 2 /A)A tk (0,q;0), (11) 

p 71 k Jo 
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where 



q /r • pq ■ q- - ////•;/// 2 

Q(p,q;E) = — log- — —2, (12) 

2p pi — pq + — m E l% 



D k (p 2 ) = [-l/a k + ^p 2 -ie\ , (13) 

and A is an ultraviolet cutoff that must be large compared to p, 1 / 1 a 2 j , and l/|a 3 |. 

Since the T-matrix elements in Eqs. (9) and (10) involve only the three linear combina- 
tions Yli a iAij(0,p; 0), the 9 coupled STM equations can be reduced to 3 coupled integral 
equations for these 3 linear combinations. If A is sufficiently large, the solutions to the inte- 
gral equations in Eqs. (11) depend only log-periodically on A with a discrete scaling factor 
Ao ~ 22.7. The dependence on the arbitrary cutoff A can be eliminated in favor of a physical 
3-body parameter, such as the Efimov parameter defined by Eq. (1). If we restrict A to 
a range that corresponds to a multiplicative factor of 22.7, then A differs from only by a 
multiplicative numerical constant. Thus, we can also simply use A as the 3-body parameter 
[36]. ( 

If A is real valued, the STM equations describe atoms that have no deep dimers. The 
rate K^ ccp for 3-body recombination into deep dimers, which can be obtained by combining 
Eqs. (5), (6), (9), and (10), must therefore be zero. For atoms that have deep dimers, the 
effects of the deep dimers can be described indirectly by using a complex-valued 3-body 
parameter. If we use the ultraviolet cutoff A as the 3-body parameter, the inclusive effects 
of deep dimers can be taken into account by analytically continuing the upper endpoint of 
the integral A in Eqs. (11) to a complex value Aexp(i?7*/so)- The path of integration in 
the variable q can be taken to run along the real axis from to A and then along the arc 
from A to A exp(ir]* / So) . This path can be deformed to run along the straight line from to 
Aexp(ii]^/s ) provided we add explicitly the contributions from any poles that are crossed 
as the contour is deformed. These poles arise from the diatom propagator D k (p 2 ), which in 
the case a k > has a pole associated with the shallow dimer at p — l/a k . For example, if 
a k > 0, the integral in Eq. (11) with a complex ultraviolet cutoff Ae ir '^ So can be written 

J dq + j dq) Q(p, q; 0)D k (3q 2 /A)Ai k (0, q\ 0) 
Alii i-Ac'W'o 

= Q(p,q k ;0)-^A ik (0,q k ;0) + J dqQ{p,q;0)D k (3q 2 /A)A ik (0,q;0), (14) 



where q k = 2/(\^3a k ). In the integral on the right side, the integration contour runs along the 
straight line path from to Ae ir, */ S °. With the ultraviolet cutoff replaced by Aexp(^/s ), 
the rate K^ ccp for 3-body recombination into deep dimers is nonzero. 

The rate constant K 3 in Eq. (9) requires the extrapolation of the solutions Aij(0,p; 0) to 
the STM equations in Eqs. (11) to p — 0. These solutions are singular as p — » 0. The singular 
terms, which are proportional to l/p 2 , 1/p, and hip, appear only in Re Aij(0, p; 0) for real 
p and can be derived by iterating the integral equations [37]. Since lmAij(0, p\ 0) must be 
extrapolated to p — 0, it is useful to transform Eqs. (11) into coupled STM equations for 
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amplitudes *Aij(0,p; 0) obtained by subtracting the singular terms from Aij(0,p; 0): 

1 S 7T 

A tJ (0,p; 0) = Aj(0,p; 0) V 1 + J2a n (l - 5 in )(l - 5 nj ) 



~ l0g X I ^S ^ 1 "^"X 1 _ M 

I n 

-g J^a m a n (l - <J im )(l - <J m „)(l - <J„j) > • (15) 

m,n J 

D. Three Equal Scattering Lengths 

We can obtain analytic results for the 3-body recombination rate in the case of three 
equal scattering lengths: a\i = 023 = = a. In this case, the 3-body recombination rates 
in Eqs. (9) and (10) depend only on the combination ^2iAij(0,p; 0)- O ne can show that the 
solutions to the STM equation for Yli^ijiPj Q'i E) are ^he same for j = 1, 2, 3 and are equal 
to the corresponding amplitude A(p, q; E) for identical bosons with scattering length a: 

J2 Aj{p, 5; E) = A{p, q; E), j = 1, 2, 3. (16) 

i 

The STM equation analogous to Eq. (11) for identical bosons is 1 

A(0,p;0) = 4 + - / dqQ(p,q;0)D(3q 2 /A)A(0,q;0). (17) 

If the low-temperature limit is taken with the number density for the identical bosons much 
smaller than the critical density for Bose-Einstein condensation, the rate constant for 3-body 
recombination of identical bosons into the shallow dimer is [8] 



K: 



shallow 
3 



512vr 2 ft 



A(0,2/(V3a);0) . (18) 

V om 

Upon making the substitution ^2 t aiAik — > ciA in Eq. (10), we see that the recombination 
rate for the three fermions into the (ij) dimer is exactly 1/3 of the recombination rate 
for identical bosons in Eq. (18). Summing over the three shallow dimers, we find that the 
expression for the recombination rate of the three fermions into shallow dimers is identical to 
the expression for the recombination rate of identical bosons into the single shallow dimer. 
Similarly, we find that the expression in Eq. (9) for the total recombination rate of the three 
fermions is identical to that for the total recombination rate of identical bosons. 

Braaten and Hammer deduced a semi-analytic expression for the rate constant for 3-body 
recombination of identical bosons into deep dimers with a large negative scattering length 
a [16]: 

if 3 dccp = 16 ^ si * h (2^) ha? 
sin [sq ln(a/a'J] + sinh 77* m 



1 The amplitude A(p,q;E) differs from the amplitude As(p,q;E) in Rcf. [8] by a multiplicative factor of 
o/(8tt). 
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where sq = 1.00624, a* = —1/(Dk*), and C and D are numerical constants. The most 
accurate values for the numerical constants are C = 29.62(1) and D = 0.6642(2) [28, 38]. 
The relation between the Efimov parameter and the ultraviolet cutoff is = 0.17609(5) A, 
modulo multiplication by an integer power of A ~ 22.7. The expression for the recombi- 
nation rate constant in Eq. (19) exhibits resonant enhancement for a near the values A^a* 
for which there is an Efimov trimer at the 3-body threshold. The line shape in Eq. (19) 
for the 3-atom loss resonance as a function of the scattering length played a key role in the 
discovery of an Efimov state for 133 Cs atoms in the lowest hyperfine state [12]. It applies 
equally well to three fermions with equal negative pair scattering lengths. 

Macek, Ovchinnikov, and Gasaneo [39] and Petrov [40] have deduced a completely an- 
alytic expression for the 3-body recombination rate constant for identical bosons with a 
large positive scattering length a in the case where there are no deep dimers. Braaten and 
Hammer generalized their result to the case where there are deep dimers by making the an- 
alytic continuation k* — > exp(ir)* / s ) in the amplitude for this process [9]. The resulting 
analytic expression for the recombination rate is [9] 



hallow 1287r 2 (47r - 3 yg) (sin 2 [s m(a/a, )] + sinh 2 ^) ha 4 

K 3 = — T27 : — 2 r w / ni ~ (a > 0) , (20) 

smh (7rs + V*) + cos^[s ln(a/a* )J m 



where a*o ~ 0.32 k" 1 . This expression exhibits minima for a near the values Aoa*o arising 
from destructive interference between two pathways for recombination. This formula applies 
equally well to three fermions with equal positive pair scattering lengths. It gives the 
recombination rate 7^| hallow i n Eq. (6), which is summed over the 3 shallow dimers. 



E. Efimov Trimers 

The transition amplitudes Aij(p,q; E) have poles in the total energy E at the energies 
E^ of the Efimov trimers. Near the pole, the amplitudes factor: 

A (n a- F) , ?MMl (01) 
Aijip, q, E) — )■ E _ E(n) ■ (21) 

The spectrum of Efimov trimers can be obtained by solving the three coupled homogeneous 
integral equations for Bj(q): 

^(p) = -y>-<^) / dqQ(p,q;E)D k (3q 2 /A-mE/h 2 )B k (q), (22) 

where the ultraviolet cutoff A must be much larger than p, \mE /Ti 2 \ 1 / 2 , and all three inverse 
scattering lengths 1/aj. The set of homogeneous STM equations in Eq. (22) are nonlinear 
eigenvalue equations for E. If A is real valued, the eigenvalues are real valued if the energy 
is below all the scattering thresholds. The Efimov trimers are therefore sharp states with 
widths. There is always a 3-atom scattering threshold at E — 0. If a^- > 0, the atom-dimer 
scattering threshold at E — —h 2 / (ma 2 j) has lower energy. 

If there are deep dimers, the Efimov trimers can decay into a deep dimer and a recoiling 
atom. Their widths can be calculated by analytically continuing the upper endpoint of the 
integral A in Eqs. (22) to a complex value Aexp(?^/s ). The complex energy eigenvalue for 
the trimer can be expressed as 

#W = -4") -irP/2. (23) 
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If is small compared to the difference between and the nearest scattering threshold, 
then and can be interpreted as the binding energy and the width of the trimer, 

respectively. If is not small, they do not have such precise interpretations. 

If ?7* <C 1, the binding energies and widths of the Efimov trimers can be calculated 
approximately by solving the STM equation with a real-valued cutoff A. The complex 
energy E^ in Eq. (23) is a function of the complex parameter Aexp(ir]^/s ) that must be 
real valued in the limit 77* — > 0. Expanding that function in powers of 77*, we find that the 
leading approximations to the binding energy and the width are 

4 B) « -£ (n) Uo> ( 24 ) 

(25) 



I>) « _2* A £(») 
T s dA 



v*=o 



The derivative with respect to A in Eq. (25) is calculated with the scattering lengths 
held fixed. The leading corrections to Eqs. (24) and (25) are suppressed by a factor of 77^. 



F. Dimer Relaxation 

If our system of fermionic atoms is a mixture of (jk) dimers and atoms of type i, another 
3-atom loss process is dimer relaxation: the inelastic scattering of the atom and the (jk) 
dimer into an atom and a dimer with a larger binding energy. If % coincides with j so there are 
only two distinct spin states, the dimer relaxation rate decreases the scattering 

length increases [41]. The atom-dimer mixture is therefore remarkably stable when the 
scattering length is very large. This stability was verified in the recent experiments with 
6 Li atoms [26, 27]. If % is distinct from both j and k so there are three distinct spin states, 
there is no Pauli suppression of the atom-dimer relaxation rate. The rate increases as the 
scattering lengths are increased. On top of that, there can also be resonant enhancement 
associated with Efimov physics. 

To be definite, we consider a mixture of (23) dimers and atoms of type 1. We denote the 
number densities of atoms and (23) dimers by ri\ and 71(23), respectively. The loss rate from 
atom-dimer relaxation can be expressed in terms of a rate constant /?i(23): 

d d 

Atom-dimer relaxation channels are also inelastic atom-dimer scattering channels. So if 
Pi(23) > 0, the atom-dimer scattering length 01(23) must have a negative imaginary part. 
These two quantities are related by the optical theorem: 

01(23) = Im(a 1(2 3)). (27) 

m 

If the scattering lengths are all large, then by dimensional analysis, /?i(23) must be 0023/772 
multiplied by a dimensionless coefficient that depends on the ratios of scattering lengths 
012/^23 and 013/023 and also on 023 re*, where «* is the Efimov parameter. The dependence 
on (I23K* is required to be log-periodic with discrete scaling factor A ~ 22.7. The universal 
predictions for dimer relaxation rates can be calculated by solving appropriate sets of coupled 
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STM equations. The dimensionless coefficient of %a2^/m in the relaxation rate constant $1(23) 
can be especially large if there is an Efimov trimer close to the l+(23) atom-dimer threshold. 
In this case, there is resonant enhancement of the dimer relaxation rate. The resulting loss 
feature is called an atom-dimer loss resonance. 

ILL 6 LI ATOMS: LOW-FIELD UNIVERSAL REGION 

In this section, we apply our formalism to the lowest three hyperfine spin states of 6 Li 
atoms in the region of low magnetic field from to 600 G. 

A. Lowest Hyperfine Spin States of 6 Li Atoms 

An example of a fermion with three spin states is 6 Li atoms in the three lowest hyperfine 
spin states. They can be labelled by their hyperfine quantum numbers \f,rrif) or by integers: 
|1) = ||,+|), |2) = ||,— |), and |3) = |§,— §)■ The pair scattering lengths ai 2 , a 23 , and 
ai3 have Feshbach resonances near 834 G, 811 G, and 690 G, respectively [25]. Beyond 
these Feshbach resonances, all three scattering lengths approach the triplet scattering length 
—2140 ao, which is large and negative. A convenient conversion constant for 6 Li atoms is 
h/m = 1.0558 x 10~ 4 cm 2 /s. 

The relevant range for ultracold atoms is the van der Waals length £ v dw = (mCs/h 2 ) 1 ^, 
which is approximately 62.5 a for 6 Li. The corresponding energy scale is the van der Waals 
energy E vd w = h 2 /(m£l dW ). The associated frequency is z/ v dw = E v dw/(2irh) 154 MHz. 
The zero-range approximation should be accurate if | a-12 1 , 1 023 1 , and 1 0,13 1 are all much larger 
than £ v dw and it should be at least qualitatively useful if |ay| > 2£ v dw- There are two regions 
of the magnetic field in which all three scattering lengths are larger than 2£ v d\y ~ 125 a : 
a low-field region 122 G < B < 485 G and a high-field region B > 608 G. These two 
universal regions are separated by a nonuniversal region in which all three scattering lengths 
go through zeros. Efimov physics in the universal regions will be characterized by values of 

and 77* that may not be the same in the two regions. In general, these parameters may 
be expected to vary slowly with the magnetic field, just like the scattering length away from 
a Feshbach resonance. In a sufficiently narrow region of magnetic field, they can be treated 
as constants. While their values could in principle be calculated from microscopic atomic 
physics, in practice they have to be determined by measurements of 3-body observables. 

In Fig. 1, the three scattering lengths ai 2 , a 2 3, and a 13 are shown as functions of the 
magnetic field in the low-energy region from to 600 G [42]. Throughout most of this 
region, the smallest scattering length is a 12 . It satisfies |a 12 | > 2 £ v dw i n the interval 
122 G < B < 485 G and achieves its largest value —290 a = —4.6 ^ v dw near 320 G. 
This interval therefore contains a universal region in which all three scattering lengths are 
negative and relatively large. The zero-range approximation should be quantitatively useful 
in the middle of this interval, but it becomes increasingly questionable as one approaches 
the edges. 

B. Three-body Recombination 

The first measurements of the 3-body recombination rate K% for the three lowest hyperfine 
spin states of 6 Li atoms were carried out by Ottenstein et al. [26] and by Huckans et al. [27]. 



11 



„ 1.0 

o 

0.5 



CO 



o 0.0 
c/) 



03 



■0.5 
■1.0 



I I I I I I 


I I I I I I I I I I I I I I I I I I 

12 

23 

13 


i i i i i 

i 
JL 


t 




— — 

i i i i i i i i i 


— / v 

^ v : 
/ 

/ - 









200 400 
B[G] 



600 



FIG. 1: (Color online) The scattering lengths in units of 10 3 ao for the three lowest hyperfine states 
of 6 Li as functions of the magnetic field B from to 600 G [42]. The two vertical dotted lines mark 
the boundaries of the region in which the absolute values of all three scattering lengths are greater 
than 2 £ v dW- 



Their results for magnetic field in the region from to 600 G are shown in Fig. 2. The results 
for K 3 in Ref. [26] from measurements of the loss rates of the three individual spin states 
have been averaged to get a single value of K 3 at each value of B. Both groups observed 
dramatic variations in K 3 with B, including a narrow loss feature near 130 G and a broader 
loss feature near 500 G. 

The narrow loss feature and the broad loss feature in the measurements of Refs. [26, 27] 
both appear near the boundaries of the region in which all three scattering lengths satisfy 
\ciij\ > 2£ v dw- The effects of the finite range of the interaction may be significant near the 
boundaries of this region. In Ref. [28], we fit the data for K 3 in this region by calculating the 
3-body recombination rate K^ cep in Eq. (9), using the magnetic field dependence of the three 
scattering lengths shown in Fig. 1 and treating A and 77* as fitting parameters. Since the 
systematic error in the normalization of K 3 was estimated to be 90% in Ref. [26] and 70% in 
Ref. [27], we only fit the shape of the data and not its normalization. A 3-parameter fit to 
the data from Ref. [26] in the region 122 G < B < 485 G with an adjustable normalization 
factor determines the 3-body parameters A = 436 clq 1 and 77* = 0.11. This value of the cutoff 
is equivalent to «;* = 76.8 Oq . These parameters and 77* determine the normalization of 
K^ eep , so the theoretical curve in Fig. 2 is absolutely normalized. The fit to the shape of the 
narrow loss feature is excellent. The normalization is also correct to within the systematic 
error in the data. However the fit predicts that K 3 should be almost constant in the middle 
of the low-field region and that there should be another narrow loss feature at its upper end 
near 500 G. These predictions are not consistent with the data in Fig. 2, which increases 
monotonically in the middle of the low-field region and has a broad loss feature near the 
upper end of this region. 

Similar results for the 3-body recombination rate of 6 Li atoms were obtained subsequently 
by two other groups [29, 30]. Naidon and Ueda used hyperspherical methods to calculate 
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FIG. 2: (Color online) Three-body recombination rate constant K3 as a function of the magnetic 
field B from to 600 G. The solid squares and dots are data points from Refs. [26] and [27], 
respectively. The curve is the absolutely normalized result for K^ cep with = 76.8 Oq 1 and 
77* = 0.11. The two vertical dotted lines mark the boundaries of the region in which the absolute 
values of all three scattering lengths are greater than 2 £ v dw- 

the recombination rate [29] . For their 3-body parameters, they used the real and imaginary 
parts of the logarithmic derivative of the hyperradial wavefunction. Schmidt, Floerchinger, 
and Wetterich used functional renormalization methods to calculate the recombination rate 
[30]. For their 3-body parameters, they used the real and imaginary parts of the detuning 
energy of a triatomic molecule. The results of both groups for K3 as a function of the 
magnetic field are similar to our results in Fig. 2. One difference is that in Refs. [29, 30] 
the recombination rate was calculated only up to an overall normalization constant that was 
determined by fitting the data. In our calculation, the absolute normalization is determined 
by the 3-body parameters and 77*. 

In Ref. [43], Wenz et al. provided an explanation for the loss feature near 500 G being 
much broader than predicted in Refs. [28-30] . They pointed out that there are deep dimers 
whose binding energies vary significantly over the low-field region. These dimers are those 
responsible for the Feshbach resonances 690 G, 811 G, and 834 G. In the high- field region, 
they are shallow dimers, but they become deep dimers in the low-field region. Their binding 
energies change over the low-field universal region by as much as a factor of 6. Wenz et al. 
assumed that contributions to the 3-body parameter 77* scale like the inverse of the binding 
energy of the deep dimer. They used the coefficient in this scaling relation as a fitting 
parameter along with the 3-body parameter equivalent to /t*. They obtained an excellent 
fit to K 3 over the entire low-field region, including the narrow loss feature, the broad loss 
feature, and the monotonic rise in between. Their assumption for the scaling of 77* with the 
binding energy of the deep dimer can be justified by an explicit calculation in a two-channel 
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FIG. 3: (Color online) Energy of the Efimov trimer as a function of the magnetic field B in the 
low-field region. The binding frequency /(2nh) (dark solid line) and the frequency T^) /{2ttK) 
associated with the width (difference between dark dashed lines) were obtained from the complex 
energy eigenvalue calculated using the parameters = 76.8 Oq 1 and ?]* = 0.11. Also shown for 
comparison are the small-77* approximation in Eqs. (24) and (25) for the binding frequency (light 
solid line) and for the frequency associated with the width (difference between light dashed lines). 
The two vertical dotted lines mark the boundaries of the region in which the absolute values of all 
three scattering lengths are greater than 2 £ v dw- 

model [44]. 

In Ref . [43] , Wenz et al. proposed a simple analytic approximation for the 3-body recom- 
bination rate in regions where all 3 scattering lengths are negative. Their approximation is 
the analytic result for equal scattering lengths in Eq. (19), with a replaced by an effective 
scattering length a m given by 

a m = - [{a\a 2 2 + a\a\ + a\a\) /3] 1/4 . (28) 

Another possible choice for an effective scattering length is the geometric mean of the three 
scattering lengths: 

a g = — |aia 2 a 3 | 1 / 3 . (29) 

We can use our universal results to test the accuracy of that approximation. We find that 
the analytic result in Eq. (19) with a replaced by the geometric mean a g in Eq. (29) is a 
significantly more accurate approximation. 

C. Efimov Trimers 

Naidon and Ueda [29] and Schmidt, Floerchinger, and Wetterich [30] calculated the bind- 
ing energy of the Efimov trimer that is responsible for the loss features in Fig. 2. Its binding 
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energy goes to near both the observed narrow loss feature near 130 G and the predicted 
narrow loss feature near 500 G. This demonstrates explicitly that both loss features arise 
from the same Efimov trimer crossing the 3-atom threshold. 

We calculate the spectrum of Efimov trimers by solving for the complex energy eigenvalues 
of the homogeneous STM equations in Eq. (22). The energies and widths of the trimers are 
obtained by expressing the complex energies in the form -E^ ] - irj° /2. We choose to 
label the Efimov trimer responsible for the narrow loss feature near 135 G by the integer 
n — 1. Our predictions for the binding energy E^p of the shallowest Efimov state are shown 
in Fig. 3 as a function of the magnetic field. The binding frequency E^p /{2irK) increases 
from at 134 G to a maximum of 9.56 MHz at 334 G and then decreases to at 494 G. Our 
maximum binding frequency agrees well with the maximum frequencies of about 10 MHz 
and about 11 MHz obtained in Refs. [29] and [30], respectively. The binding frequency 
E^ 1 /{2ith) associated with the next deeper Efimov trimer is predicted to change gradually 
from 12.1 GHz at 134 G to 12.5 GHz at 334 G and then to 12.1 GHz at 494 G. Since this 
binding energy is much larger than the van der Waals frequency z/ vdW = 154 MHz, this 
Efimov trimer and all the deeper ones predicted by the STM equations are artifacts of the 
zero-range approximation. Our value of /t* can be interpreted as the binding wavenumber 
{mEy ^ /h 2 ) 1 / 2 of the fictitious Efimov trimer labelled by n = —2. It corresponds to the 
choice n* = —2 in Eq. (1). 

We note that the finite range corrections to the binding energies of the deeper trimers 
can be substantial if their energy becomes comparable to the van der Waals energy even if 
all scattering lengths are large. The leading corrections to the binding energy for the deep 
trimers are of order £ vdW (mE^ /h 2 ) 1 / 2 . For the trimer labeled n — 1, these corrections are 
25% for the largest value of the binding energy. 

In Ref. [29], Naidon and Ueda also calculated the width of the Efimov trimer. In Fig. 3, 
our results for the width are illustrated by plotting the frequencies associated with the 

energies E^ ±T^/2 as functions of the magnetic field. The frequency /(27rh) increases 
from 2.73 MHz at 134 G to a maximum of 7.14 MHz at 332 G and then decreases to 2.56 MHz 
at 494 G. Our maximum width is more than twice as large as the maximum width of about 
(3 MHz)x27r/i obtained in Ref. [29]. Since our result for is always comparable to or 
larger than E^ , the interpretations of E^ and as the binding energy and width of the 
trimer should be viewed with caution. 

Since i]* = 0.11 is relatively small, we can also use the small- 1]* approximations to the 
binding energy and the width given in Eqs. (25) and (24). These approximations, which 
are calculated using = 76.8 (Xq ; ctrc shown in Fig. 3 for comparison. The maximum 
binding frequency E!p / (2k%) is 10.1 MHz at 334 G, which is larger than the result from the 

complex energy by about 6%. The frequency /{2nh) increases from 2.76 MHz at 134 G 
to a maximum of 7.23 MHz at 332 G and then decreases to 2.59 MHz at 494 G. These values 
are larger than the results from the complex energy by only about 1%. 

IV. 6 LI ATOMS: HIGH-FIELD UNIVERSAL REGION 

In this section, we apply our formalism to the lowest three hyperfine spin states of 6 Li 
atoms in the region of high magnetic field from 600 to 1200 G. In Fig. 4, the three scat- 
tering lengths ai2, 023, and a i3 are shown as functions of the magnetic field [42]. This 
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FIG. 4: (Color online) The scattering lengths in units of 10 3 ao for the three lowest hyperfine states 
of 6 Li as functions of the magnetic field B from 600 G to 1200 G [42]. The three vertical lines 
mark the positions of the Feshbach resonances. 



region includes the Feshbach resonances in a 12 , 023, and a 13 near 834 G, 811 G, and 690 G, 
respectively. Beyond these Feshbach resonances, all three scattering lengths approach the 
triplet scattering length —2140 a . For B > 637 G, the absolute values of all three scat- 
tering lengths are larger than 2140 a® ~ 34 £ v dw- An estimate of the lower boundary of 
this universal region is 608 G, where the smallest scattering length is 013 = 125 a ~ 2£ v dw- 
The zero-range approximation should be very accurate throughout most of this high-field 
region. The physics in this universal region is rich, with the Feshbach resonances marking 
the boundaries between regions in which 3, 2, 1, or of the three scattering lengths are 
positive. 



A. Measurements of Three-body Recombination 

In Ref . [27] , Huckans et al. presented measurements of the 3-body recombination rate at 
6 values of the magnetic field in the range 600 G to 1000 G. The low-temperature limit of the 
recombination rate could not be measured, because the temperature was not low enough 
and because of heating associated with the Feshbach resonance. In Ref. [28], we showed 
that a naive fit of the last two data points at 894 G and 953 G using our zero-temperature 
calculations indicated an Efimov resonance near 1200 G. We suggested that it might be 
worthwhile to search the high-field region for an Efimov resonance. 

A narrow 3-atom loss resonance in the high-field region near 895 G was recently discovered 
by Williams et al. [31] and by Jochim and coworkers [32]. Williams et al. measured the 3- 
body recombination rate at magnetic fields from 842 G up to 1500 G at temperatures lower 
than 180 nK and from 834 G up to 955 G at temperatures lower than 30 nK. Their data 
is shown in Fig. 5. By fitting their measurements using our formalism, they obtained the 
3-body parameters = (6.9 ± 0.2) x 10~ 3 a^ 1 and 77* = 0.016±g;g$g [31]. Since «, is 
only defined up to multiplication by integer powers of Aq ~ 22.7, an equivalent value is 
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FIG. 5: (Color online) Three-body recombination rate constant K3 as a function of the magnetic 
field B from 600 G to 1200 G. The solid dots and triangles are data points from Ref. [27] at 
temperatures less than 180 11K and less than 30 nK, respectively. The curve between 834 G and 
1200 G is a fit to that data. In the region from 600 G and 834 G, the curves are our predictions for 
the total 3-body recombination rate (thick line) and the contribution from recombination into deep 
dimers (thin line) for = 80.7 a^ 1 and 77* = 0.016. The three vertical lines mark the positions of 
the Feshbach resonances. 

= 80.7 ± 2.3 <2q . Their fit for the central values of and 77* is shown in Fig. 5 as 
the thick solid line between 834 G and 1200 G. The fitted value for the position of the 
3-atom loss resonance is 895^ G. The peak value of the 3-body recombination rate at zero 
temperature is predicted to be x 10~ 17 cm 6 /s. 

It is interesting to compare the values of the 3-body parameters and 77* in the high- 
field region with those in the low-field region. The values in the low-field region obtained 
by fitting the measurements of the narrow 3-body recombination loss feature by Ottenstein 
et al. [26] were 76.8 a^ 1 and 77* ~ 0.11. It is difficult to quantify the errors in these 

parameters, because the narrow loss feature is at the edge of the universal region where 
range corrections may be significant. The values of k* in the high- field region and the low- 
field region are consistent within errors. This could be just a coincidence, but it suggests 
that the phase in the 3-body wavefunction at short distances that controls Efimov physics 
is insensitive to the Feshbach resonances that change the scattering lengths. The value 
of 77* in the high-field region is about an order-of-magnitude smaller than in the low-field 
region. This demonstrates that the 3-body parameters need not be equal in the two universal 
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FIG. 6: (Color online) Energies of the Efimov trimers as functions of the magnetic field B in 
the high-field region. The solid curves are the predicted binding frequencies E^/(2irh) for = 
80.7 1 and rj* = 0.016. The dashed curves are the upper and lower error bars obtained by varying 
K*. The curves labelled 12, 23, and 13 are the atom-dimer thresholds. The dots indicate the points 
where the trimers disappear through the l+(23) atom-dimer threshold. The horizontal dotted line 
is the van der Waals frequency 154 MHz. 

regions. Wenz et al. proposed a mechanism for variations in 77*, namely significant changes 
in the binding energies of deep dimers with the magnetic field [43]. This mechanism might 
also explain the order-of-magnitude difference in 77* between the two regions. 

Using the values of and 77* determined by O'Hara et al. [31], we can calculate the 
universal predictions for other aspects of Efimov physics for 6 Li atoms in the high-field 
universal region. The error bars on the values of and 77* can be used to give error bars 
on the predictions. In the next subsections, we give universal predictions for the binding 
energies and the widths of the Efimov trimers and for the 3-body recombination rate. 

B. Efimov Trimers 

We calculate the binding energies and the widths of the Efimov trimers by solving for 
the complex energy eigenvalues of the homogeneous STM equations in Eq. (22). We choose 
to label the Efimov trimer responsible for the narrow loss feature near 895 G by the integer 
n — 1. In Fig. 6, the predicted binding energies E^p and E^p of the two shallowest Efimov 
trimers are shown as functions of the magnetic field. The Efimov trimer responsible for the 
narrow loss feature has a binding energy that vanishes at 895^ G. As the magnetic field 
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FIG. 7: Widths of the Efimov trimers as functions of the magnetic field B in the high-field region. 
The solid curves are the frequencies Tjf* /(2nh) associated with the widths for rc* = 80.7 a^ 1 and 
77* = 0.016. The dashed curves are the upper and lower error bars obtained by varying ry*. 

decreases, increases monotonically until the critical value = 672 ± 2 G where the 
Efimov trimer disappears through the l + (23) atom-dimer threshold. When it disappears, its 
binding energy relative to the 3-atom threshold is 871+gg kHz, which is also the binding 
energy h 2 / (rna 2 3 ) of the (23) dimer. The next deeper Efimov trimer has a binding frequency 
EP /(2nh) that increases monotonically from 26.4^ MHz at 895 G to 34.9^11 MHz at B*. 
Our solutions to the STM equations predict that it crosses the l+(23) atom-dimer threshold 
at 597 G, where its binding frequency relative to the 3-atom threshold is 203 MHz. This 
frequency is larger than the van der Waals frequency 154 MHz. Moreover the smallest 
scattering length at this point is only 54 ao, which is smaller than the van der Waals length 
65 ao- The zero-range predictions are not expected to be accurate for such large energies 
and for such small scattering lengths. All the deeper Efimov trimers have binding energies 
that are much larger than the van der Waals energy and are therefore artifacts of the 
zero-range approximation. Our value of /t* can be interpreted as the binding wavenumber 
(mEif^/h 2 ) 1 / 2 of the fictitious Efimov trimer labelled by n = —2. It corresponds to the 
choice n* = —2 in Eq. (1). 

As discussed in Sec. Ill C, the leading range corrections to the binding energy for the 
deep trimers are estimated to be of order l vdW (mE^ /Ti 2 ) 1 / 2 . For the trimer labeled n = 1, 
these corrections are 8% where it crosses the l+(23) atom-dimer threshold and become 
significantly smaller at larger values of the magnetic field. For the trimer labelled n = 
these corrections are about 40% at B* and larger values of the magnetic field. 

In Fig. 7, the predicted widths of the Efimov trimers are shown as functions of the 
magnetic field. The frequency /{2'kK) associated with the width of the shallower trimer 
increases from 0.706 kHz at 895 G to a maximum of 7.98 kHz at 692 G and then decreases to 
0.940 kHz at The frequency /(27rh) associated with the width of the deeper trimer 
increases from 1.74 MHz at 895 G to 2.00 MHz at B*. 
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FIG. 8: (Color online) Three-body recombination rate constant K3 as a function of the magnetic 
field B from 650 G to 850 G. The curves are our predictions for the total 3-body recombination 
rate (thick solid line), the contribution from recombination into deep dimers (thin solid line), 
and the contributions from recombination into the (12) dimer (short-dashed line), the (23) dimer 
(dash-dotted line), and the (13) dimer (long-dashed line) for = 80.7 clq 1 and 77* = 0.016. 

C. Predictions for Three-body Recombination 

Our predictions for the 3-body recombination rate in the region 600 G to 834 G, where 
one or more of the scattering lengths are positive, are shown in Fig. 5. The upper curve is 
the total recombination rate K 3 , while the lower one is the contribution i^3 Gep from deep 
dimers. Besides the Efimov resonance at 895 G, the only other peaks are at the Feshbach 
resonances at 690 G, 811 G, and 834 G. They arise simply from the scaling of K3 as a 4 up to 
logarithms, where a is some appropriate mean of the scattering lengths a\2, 0,23, and 013. The 
total recombination rate K$ is predicted to have three local minima. There is a broad local 
minimum at 759 ± 1 G, where the minimum recombination rate is (2.51q; 5 ) x 10~ 22 cm 6 /s. 
There is a narrow local minimum at 829 ± 1 G, where the minimum recombination rate is 
(2.2±H) x 10~ 18 cm 6 /s. Finally there is a local minimum at 861 ± 2 G between the last 
Feshbach resonance and the 3-atom Efimov resonance, where the minimum recombination 
rate is (5.8±|;g) x 10" 19 cm 6 /s. 

In Fig. 8, we show the predictions for the 3-body recombination rate in the region from 
650 G to 850 G in more detail. In addition to K3 and K^ eep , we show the contributions 
from recombination into the (12), (23), and (13) dimers. This figure reveals that the local 
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minima at 759 G and 829 G are associated with interference in the recombination into 
shallow dimers. The narrow minimum at 829 G arises from the combination of a rapidly 
increasing rate into deep dimers and a rapidly decreasing rate into (12) dimers, which has its 
minimum at 830 G. The broad minimum at 759 G arises from minima in the rates into (12) 
and (23) dimers at 757 G and 765 G, respectively. There are also minima in the rates into 
the (13) dimer at 672 G and the (23) dimer at 600 G, but their effects are not visible in the 
total recombination rate. We have verified that these minima in the recombination rates into 
individual shallow dimers arise from interference effects by showing that they become zeroes 
as r?* is decreased to 0. Thus these minima are the result of destructive interference between 
two recombination pathways. In the case of three equal positive scattering lengths, a similar 
interference effect is evident in the analytic expression for the 3-body recombination rate 
into shallow dimers in Eq. (20). 

D. Atom-Dimer Resonance 

An atom-dimer loss resonance can appear at a value of the scattering length for which 
an Efimov trimer crosses the atom-dimer threshold. From the Efimov trimer spectrum 
in Fig. 6, one can see that atom-dimer resonances are predicted at the two values of the 
magnetic field where the Efimov trimers cross the l + (23) atom-dimer threshold. The atom- 
dimer resonance associated with the shallower of the two Efimov trimers is predicted to 
occur at B* = 672 ±2 G. The universal predictions for the dimer relaxation rate near this 
resonance can be calculated by solving appropriate sets of coupled STM equations. The 
atom-dimer resonance associated with the deeper of the two Efimov trimers is predicted 
to occur at B'^ = 597 G. This resonance occurs slightly outside the universal region, so 
universal predictions for the position of the resonance and for the dimer relaxation rate are 
not expected to be accurate. 

If we restrict our attention to magnetic fields very near the atom-dimer resonance at B*, 
we can get an approximation to the universal predictions for the dimer relaxation rate /?i(23) 
without actually solving the STM equations. We take advantage of the fact that the atom- 
dimer scattering length au2z) diverges at B*. The 3- atom problem in this region therefore 
reduces to a universal 2-body problem for the atom and the (12) dimer. The universal 
properties are determined by the large scattering length ai^3)- For example, in the region 
B > B* but very close to B*, the binding energy of the Efimov trimer is well approximated 
by the sum of the binding energy of the (23) dimer and a universal term determined by the 
atom-dimer scattering length: 



This should be a good approximation as long as the second term is much smaller than the 
first term. In the region near B*, the atom-dimer scattering length can be approximated by 
an expression that has the same form as the analytic result for identical bosons but with 
different numerical coefficients: 



where a* = 023 (-B*). The coefficients C\ and C2 can in principle depend on ratios of the 
scattering lengths, but 013 is an order of magnitude larger so it essentially decouples and a±3 




(30) 



01(23) ~ (Cicot[soln(a 23 /a*) + i?7*] +C 2 )a 23 , 



(31) 
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FIG. 9: (Color online) Dimer relaxation rate constant /?i(23) for (23) dimers and atoms of type 1 
as a function of the magnetic field B. The magnetic field is measured relative to the position 
of the atom-dimer resonance. The curves are predictions using the approximation in Eq. (32) for 
three values of 77*: 0.006, 0.016, and 0.022. 



and 023 are both increasing with B so their ratio does not change rapidly near B*. Thus 
we can treat C% and C2 as numerical constants. These constants can be determined by 
using Eq. (30) to extract 01(23) from the results for the real part of the trimer binding energy 
shown in Fig. 6 and then fitting those results to the expression for the atom-dimer scattering 
length in Eq. (31). The resulting constants are C% = 0.67 and C2 = 0.65. 

With the constant C\ in hand, we can obtain an approximation for the dimer relaxation 
rate constant /3j(23) f° r B near B* simply by inserting the approximation for 01(23) m Eq. (31) 
into the optical theorem relation in Eq. (27). The resulting expression is 

9.44Cisinh(27fc) ha 23 

Pi(23) ~ - or — ; — ; ; — r; rri • (32) 

sin z [soln(a23/a*)J + smh 77* m 

If = 672 G, the value of a* is 835 a . The only dependence on the magnetic field is 
through the .B-dependence of 023. This approximation for /3i(23) is shown as a function of 
B — B* in Fig. 9. The maximum value of /?i(23) at the peak of the resonance is 4 x 10~ 7 cm 3 /s 
for 77* = 0.016. 

Since the atom-dimer resonance associated with the deeper of the two Efimov trimers 
occurs slightly outside the universal region, universal predictions for the dimer relaxation 
rate are not expected to be accurate. They might however be useful for making order-of- 
magnitude estimates. If we apply the analysis described above to the binding energy of the 
deeper Efimov trimer, we get the constants C\ = 0.7 and C2 = 0.3. The value of a* is 



22 



approximately 54 do- If we use Eq. (32) to estimate the dimer relaxation rate, we obtain a 
maximum value of ^1(23) at the peak of the resonance of 2.4 x 1CT 8 cm 3 /s for 77* = 0.016. 

E. Many-body physics 

One of the most important motivations for understanding few-body physics is that there 
are aspects of many-body physics that are controlled by few-body observables. Our universal 
results for the few-body physics of 6 Li atoms with 3 spin states provides useful information 
about many-body systems of 6 Li atoms with sufficiently low temperature and sufficiently low 
number densities. If the atoms are in thermal equilibrium at temperature T, an important 
scale is the thermal length \ T = (2irnik B T /h 2 )~ l l 2 . If the atoms of type % have number 
density rij, another important scale is the Fermi wavenumber kp% = (Sn 2 ^) 1 ^ 3 . In our 
calculation of the 3-body recombination rate, we set the wavenumbers of the incoming 
atoms to 0. In order for our result to apply quantitatively to all atoms in the system, it 
is necessary that A T X , kpi, and fcp 3 all be small compared to the wavenumber scale 
(l/a? + l/a 2 , + l/al) 1 / 2 set by the interactions. Even if this condition is not well satisfied, 
our result can be used to obtain order-of-magnitude estimates of the relevant time scales. 

There are several relevant time scales that we can extract from our few-body calculations. 
If the system contains Efimov trimers, an important time scale is the lifetime of the 

Efimov trimer, where is its width. The universal predictions for the widths of the Efimov 
trimers are shown in Fig. 7. If the system can be approximated by a gas of individual low- 
energy atoms, the time scale for disappearance of a significant fraction of the atoms is set 
by the 3-body recombination rate constant K%. From the rate equation in Eq. (4), we see 
that the time scale for loss of a significant fraction of the atoms of type i is {K^rijUk)^ 1 , 
where j and k are the two complimentary spin states. The universal predictions for the 
3-body recombination rate are shown in Figs. 5 and 8. Finally if the system contains 
low-energy dimers and low-energy atoms in the complimentary spin state, the time scale for 
disappearance of a significant fraction of the atoms or dimers is set by the appropriate dimer 
relaxation rate constant. From the rate equation in Eq. (26), we see that the time scales 
for loss of significant fractions of (23) dimers and of atoms of type 1 are (/3i(23)^i) _1 and 
(/3i(23)^(23)) -1 , respectively. The only universal information we have about dimer relaxation 
rates is the approximation for ^1(23) in Eq. (32), which should be accurate for magnetic fields 
within about 10 G of m 672 G. For magnetic fields further from this critical value but 
still within the region where 023 is positive, one may be able to use the extrapolation of this 
expression as an estimate for the dimer relaxation rate. 

If the system contains a Bose-Einstein condensate of low-energy dimers as well as low- 
energy atoms in the complimentary spin state, another important few-body observable is 
the atom-dimer scattering length for that particular atom and dimer. The real part of the 
atom-dimer scattering length determines the mean-field energy of the atom in the dimer 
condensate. The mean-field energy of an atom of type 1 in the (23) dimer condensate with 
sufficiently low number density 71(23) is (37rh 2 /m) Re(ai(23))7i(23)- The approximation for 
01(23) in Eq. (31) should be accurate for magnetic fields within about 10 G of B* ss 672 G. 

We will illustrate the relevance of our universal result to many-body physics by applying 
them to two specific values of the scattering length that are interesting from a symmetry 
perspective. If the pair scattering lengths and are equal, atoms of types j and k are 
related by an SU (2) symmetry. If all three scattering lengths are equal, the three spin states 
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are related by an SU(3) symmetry. There is an SU(2) symmetry point at 731 G, where 
the scattering lengths are a\2 = 0,23 ~ +2500 ao and 013 ~ —7100 ao- An SU (3) symmetry 
point can be approached by going to very high magnetic fields, where all three scattering 
lengths approach the spin-triplet scattering length —2140 a [25]. 

A quantum degenerate Fermi gas of 6 Li atoms with approximately equal populations of 
the three lowest hyperfine spin states at 1500 G has been realized by Williams et al. [31]. 
The scattering lengths are Oi 2 ~ —2460 a , a 2 3 ~ —2360 a , and a 13 ps —2240 a , so there 
is an approximate SU (3) symmetry. One candidate for a metastable ground state is a state 
that can be approximated by filled Fermi spheres for all three spin states, with Cooper 
pairing that breaks the SU(3) symmetry down to a U(l) subgroup. The lifetime for such a 
system is determined by the 3-body recombination rate K3, which is predicted to be about 
8 x 10~ 22 cm 6 /s at 1500 G. Another candidate for a metastable ground state is a filled Fermi 
sphere of Efimov trimers. The binding energy of the Efimov trimer is predicted to be 
about 25 MHzx(2tto) at 1500 G. Its width rj? is predicted to be about 1.7 MHzx(2tto), 

which corresponds to a lifetime h/T^ of about 9 x 10~ 8 s. 

We now consider a many-body system at 731 G, where there is an SU(2) symmetry 
relating the atoms of types 1 and 3. For a state that can be approximated by filled Fermi 
spheres for all three spin states, the lifetime is determined by the the 3-body recombination 
rate K 3 , which is predicted to be about 7 x 10 -22 cm 6 /s at 731 G. A better candidate 
for a metastable ground state is a state containing a Bose-Einstein condensate of dimers, 
which breaks the SU(2) symmetry down to a U(l) subgroup, and a filled Fermi sphere of 
the complimentary atoms. The lifetime of the state is determined by the dimer relaxation 
rates /?i(23) and /3 3 (i2), which are equal by the 577(2) symmetry. We can estimate /?i(23) 
by extrapolating the expression for the dimer relaxation rate /?i(23) in Eq. (32) to 731 G, 
which gives 3 x 10 -10 cm 3 /s. We can also estimate the mean-field energy of an atom of 
type 1 in a (23) dimer condensate with number density 71(23) by extrapolating the expression 
for the atom-dimer scattering length 01(23) m Eq. (31) to 731 G. The resulting estimate is 
2 x 10~ 9 Hz cm 3 x (27i7i 71(23)). The positive sign of the real part of 01(23) implies that the 
atoms of type 1 are repelled by the (23) dimer condensate. A state in which the dimer 
condensate and the atoms are spatially separated is therefore energetically favored over a 
homogeneous state. A final candidate for a metastable ground state is a filled Fermi sphere 
of Efimov trimers. The binding energy E^p of the shallower Efimov trimer is predicted to 
be about 190 kHzx(2?ro) at 731 G. Its width r^ 1} is predicted to be about 6 kHzx(2vro) at 
731 G, which corresponds to a lifetime h/F^p of about 3 x 10 -5 s. 

V. SUMMARY 

Systems consisting of 6 Li atoms with 3 spin states provide a rich playground for the 
interplay between few-body physics and many-body physics. The experimental study of 
many-body physics is only possible if the loss rates of atoms from few-body processes are 
sufficiently low. Measurements of the position and width of a single Efimov loss feature can 
be used to determine the 3-body parameters and 77*. Calculations in the zero-range limit 
can then be used to predict few-body reaction rates in the entire universal region. 

In the low-field region for 6 Li atoms, the universal predictions were only qualitatively 
successful. A fit to the measurements of the 3-body recombination rate by Ottenstein et 
al. gives the 3-body parameters /c* 77 Oq 1 and 77* ps 0.11. With these parameters, the 
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universal results for the 3-body recombination rate as a function of the magnetic field give 
a good fit to the narrow loss feature near 130 G but do not agree well with measurements 
in the upper half of the low- field region [28-30] . A reasonable explanation was proposed by 
Wenz et al. [43] : 77* is particularly sensitive to the binding energies of the shallowest of the 
deep dimers and there are 6 Li dimers whose binding energies change dramatically across the 
low-field region. Assuming the scaling behavior 77* ~ E^ cp , they obtained a good fit to the 
recombination rate in the low-field region. 

In the high-field region for 6 Li atoms, the scattering lengths are much larger so the 
universal predictions should be much more accurate. A narrow 3-atom loss feature near 
895 G was discovered by Williams et al. [31] and by Jochim and coworkers [32]. By fitting 
their measurements of the 3-body recombination rate, Williams et al. determined the 3-body 
parameters associated with Efimov physics to be /t* = 80.7 ± 2.3 1 and 77* = O.OIG^q^o- 
We used those parameters to calculate the universal predictions for the binding energies and 
widths of the Efimov trimers shown in Fig. 6 and 7. The Efimov trimer responsible for the 
narrow loss feature is predicted to disappear through the l + (23) atom-dimer threshold at 
672 ±2 G, producing a spectacular atom-dimer loss resonance. There is also a deeper Efimov 
trimer whose binding frequency and width in the universal region are approximately 30 MHz 
and 2 MHzx (2irh), respectively. This trimer is also predicted to disappear through the 
l+(23) atom-dimer threshold, but this happens outside the universal region. We also used 
the 3-body parameters determined by Williams et al. to calculate the universal predictions 
for the 3-body recombination rate, which are shown in Figs. 5 and 8. Local minima in K 3 
are predicted at 759 ± 1 G, 829 ± 1 G, and 861 ± 1 G. Finally an approximate calculation of 
the dimer relaxation rate in the region of the atom-dimer resonance is presented in Fig. 9. 
We look forward to the experimental verification of these predictions. 

In order to understand the behavior of atom-dimer mixtures at low temperatures, it 
would be useful to have universal predictions for other 3-atom observables. They include 
the atom-dimer scattering lengths, which are in general complex. The real part of an atom- 
dimer scattering length determines the mean-field shifts of the atom in the dimer condensate. 
Its imaginary part is proportional to the dimer relaxation rate. For the shallow dimer with 
the largest binding energy, which is the (23) dimer for B < 730 G and the (12) dimer for 
730 G < B < 834 G, the only relaxation channels are into deep dimers. For shallow dimers 
with smaller binding energy, there are also relaxation channels into other shallow dimers. 
It would be especially useful to have definitive universal predictions for the relaxation rate 
constant ^1(23) near the predicted atom-dimer threshold at 672 ± 2 G. It would improve 
upon the approximation illustrated in Fig. 9 by taking into account the 5-dependence of 
the ratios of scattering lengths. 

There have not yet been any direct observations of Efimov trimers in ultracold atoms. 
They have only been observed indirectly through the resonant enhancement of 3-body re- 
combination and through the resonant enhancement of atom-dimer relaxation provided by 
virtual Efimov trimers. The direct production of Efimov trimers would be another milestone 
in the study of Efimov physics in ultracold atoms. Of course, once they are produced, they 
would decay quickly. In the high-field universal region for 6 Li atoms, the deeper Efimov 
trimer is predicted to have a very short lifetime of about 10~ 7 s. The shallower Efimov 
trimer is predicted to have a lifetime of 10" 4 s to 10~ 5 s. Our universal predictions for the 
binding energies of these Efimov trimers should be useful in devising experimental strategies 
for producing them. 
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